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ABSTRACT 

The MeV radiation front of gamma-ray bursts creates copious pairs as it 
propagates through an ambient medium. The created pairs enrich the leptonic 
component of the medium by a large factor at distances R < i?ioad ~ 10^® cm 
from the burst center. The following blast wave sweeps up the pair-rich medium 
and then emits the observed afterglow radiation. We find that the afterglow 
has a “memory” of loading outside Rioad- The remain in the swept-up 
material and slowly cool down by emitting synchrotron radiation. They are likely 
to dominate the blast-wave emission in IR, optical, and UV bands during the 
hrst minutes of the observed afterglow. The expected radiation is described 
by a simple formula, which is derived analytically and checked by numerical 
integration of synchrotron emission over the blast material; a suitable Lagrangian 
formalism is developed for such calculations. The main signature of radiation 
is its flat (“white”) spectrum in a broad range of frequencies from IR to UV and 
possibly soft X-rays. This radiation can be detected by Swift satellite, which 
would enable new observational tests for the explosion physics. 

Subject headings: cosmology: miscellaneous — gamma rays: bursts — radiation 
mechanisms: nonthermal —- shock waves 


1. Introduction 

Cosmological gamma-ray bursts (GRBs) are produced by powerful explosions in distant 
galaxies. It is not yet clear how the explosion is triggered, however its basic phenomenological 
picture has been established: an ultra-relativistic shell (“fireball”) is ejected by a compact 


^Also at Astro-Space Center, Lebedev Physical Institute, Profsojuznaja 84/32, Moscow 117810, Russia 
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central engine. The expanding shell emits the burst of 7 -rays, then sweeps up an ambient 
medium, and decelerates, producing the observed afterglow radiation. The afterglow is 
explained as synchrotron emission of nonthermal electrons in the relativistic blast wave (see 
Piran 2004 for a recent review). 

Thompson & Madau (2000) pointed out that an external medium must be e^-loaded 
and preaccelerated by the leading 7 -ray front (prompt GRB radiation), which should affect 
the ensuing shock wave. The effect can be thought of as a result of 7 -ray transfer through the 
optically thin medium (Beloborodov 2002, hereafter B02) which involves a runaway of pair 
creation and bulk acceleration. A tiny fraction of GRB radiation participates in this transfer, 
however, it impacts dramatically the circumburst medium. The transfer problem was solved 
in B02, and the number of loaded and the Lorentz factor of the medium behind the 7 -ray 
front were calculated. The loaded pairs were found to dominate the ambient medium at 
radii R < i?ioad ~ 1-6 x 10^®(i?..y/10^^erg) cm, where is the isotropic ( 47 r) equivalent of 
the GRB energy. 

The loading sets the stage for the immediately following shock wave driven by the 
GRB ejecta. The shock heats and accelerates the particles of the medium, and a nonthermal 
population is expected to form behind the shock front, which produces synchrotron radi¬ 
ation. The goal of the present paper is to calculate emission from the pair-loaded postshock 
plasma in the expanding blast wave. 

This difficult problem was previously approached in a few works. B02 evaluated the 
shock parameters at i? < Rioad and found that the GRB afterglow should start with a brief 
and bright optical signal, however did not calculate the expected light curve or spectrum 
of radiation. Then Li, Dai, & Lu (2003) calculated the light curve considering the blast 
wave as a single shell with an averaged density and a common electron spectrum. This 
is not a good approximation as discussed in detail below. In particular, at i? > Rioad, only a 
small fraction of the blast-wave material is dominated by pairs, and its emission dramatically 
differs from the rest of swept-up material. 

Most recently, Kumar & Panaitescu (2004) studied e^-loaded blast waves. They focused 
on GRBs where the 7 -ray front only partially overtakes the shock wave at Rioad, so that only 
a leading part of the front creates ahead of the shock, which reduces the number of 
shocked e^. A dramatic effect appears in that situation, which was neglected in Kumar & 
Panaitescu (2004): since the postshock still overlap with the 7 -ray front, they are exposed 
to 0.1 — 1 MeV photons (keV in the plasma frame) and quickly cooled by inverse Gompton 
scattering. Most of the energy is then emitted by upscattering 0.1 — 1 MeV photons to 
GeV-TeV band, and their optical synchrotron emission is suppressed (Beloborodov 2005). 
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Whether the 7 -ray front still overlaps with the blast wave at a radius of interest depends 
on the front thickness (proportional to the duration of the prompt GRB) and the blast-wave 
Lorentz factor Fd So, two qualitatively different regimes of early afterglows are possible: 

1. — “Long-burst” (or “thick-shell”) regime where the early blast wave overlaps with the 
prompt GRB radiation. Then a strong GeV-TeV flash should be produced, and the early 
optical emission is suppressed. 

2 . — “Short-burst” (or “thin-shell”) regime where the prompt 7 -rays early overtake the 
external shock wave, and creation and Gompton cooling take place ahead of the shock. 
Then the postshock plasma is not exposed to the prompt radiation and not Gompton cooled. 
A bright optical emission can be expected in this situation. 


In the present paper, we focus on the short-burst regime. GRBs that satisfy this condi¬ 
tion have short durations th and/or modest Lorentz factors of the blast wave T (B02), 


tb < 10 



( Race \ 

lO^^cm / 


( 1 +^) 


S, 


( 1 ) 


where Race = 5“^RRioad is a characteristic radius where most of the optical-emitting pairs 
are created, and 2 ; is the cosmological redshift of the burst. We develop a suitable Lagrangian 
formalism that describes the synchrotron emission of pair-loaded blast waves and calculate 
the expected light curve and spectrum of the early afterglow. We hnd two distinct emission 
components produced at R > Rioad^ (1) the relict radiation component dominating the 
early emission from IR to soft X-rays, and (2) the recently shocked pair-free component that 
gives a standard afterglow emission initially peaking in X-rays and later evolving to softer 
bands. 


In § 2 we briefly describe the pair creation process and the formation of blast wave in the 
pair-loaded medium behind the 7 -ray front (details are found in B02). In § 3 we formulate 
the emission problem for e^-loaded blast waves and develop their Lagrangian description. 
Numerically calculated examples of e^-loaded afterglows are given in § 4. 

In § 5 we show that the component of afterglow emission is described by a simple 
formula and practically independent of the details of the shock wave physics. The case of a 
uniform ambient medium is elaborated in § 6 . In the present paper, we focus on explosions 
in media of modest density no = 0.1 — 10^ cm“^. Explosions in a high-density wind from 


^The 7 -ray front can be emitted when the explosion has a small radius, well before the blast wave forms, 
however, they may still overlap at large radii because the ejecta expands with almost speed of light. The 
7 -ray front is faster by a small Sv = cj2T'^ ~ 10“®c and completely overtakes the relativistic blast wave at 
time R/c ^ A/Sv. 
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a Wolf-Rayet progenitor will be considered in detail elsewhere. In that case, the e^-loading 
has a much stronger effect on the afterglow emission (B02). 

There are two shock fronts in a blast wave: forward and reverse. The reverse shock 
emits one more component of the early afterglow, which depends on the nature of the GRB 
ejecta. Differences between emissions from the reverse shock and the pair-loaded forward 
shock are discussed in § 7. The differences are signihcant and may allow one to distinguish 
observationally the two emission mechanisms when the Swift satellite provides the early 
afterglow data. 


2. Pair loading by the y-ray front 

A medium overtaken by a front of collimated y-rays is inevitably e^-loaded. This 
happens because some y-rays Compton scatter off the medium and get absorbed by the 
primary collimated radiation via reaction y -|- y —> e"'' + e“. 

The medium is optically thin, so only a tiny fraction of the GRB radiation front scatters 
and turns into e^, however, the number of created per ambient electron can be very 
large, n±/nQ ^ 1. The column density of an expanding photon front scales with radius as 
R~^] hence the number of scattered photons per ambient electron is decreasing as R~^ and 
n±/no 3> 1 at small R. The created do more scattering, which leads to an exponential 
runaway of pair creation. There is a sharp boundary Rioad between the exponentially loaded 
e± ];egion and the outer pair-free region. The pair creation and scattering of GRB radiation 
inside Rioad is accompanied by momentum deposition and the medium is accelerated radially 
away from the center. 

The pair loading factor n±/no and acceleration of an optically thin medium in the y-ray 
front do not depend on uq and can be calculated starting with just one ambient electron and 
one ambient proton to which the electron is coupled. The e and p components of the medium 
move always together to maintain neutrality of the plasma. Any momentum communicated 
to e is immediately shared with p, so that the effective mass of e is rUe + nip. As the y- 
ray front passes through, the electron component is enriched by a number of additional 
(then the mass per lepton is reduced) and altogether they acquire a Lorentz factor y. This 
transformation is quick: it takes time ~ (A/cy^) in the fixed lab frame, during which the 
y-ray front overtakes a given ambient electron.^ The thickness of the y-ray front A is related 
to the observed duration of the prompt GRB Rhy R = {1 + z)/S./c. 


^Velocity difference between the radiation front and the accelerated medium is (1 — 0)c « c/2y‘^. 
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A schematic explosion picture is shown in Figure 1. The radiation front leads the forward 
shock by a small distance^ 


R 

4f2 


< R, 


( 2 ) 


and changes the ambient medium just before it is shocked. The medium ahead of the blast 
wave (but already overtaken by the y-ray front) is described by the lepton number per 
ambient proton Z and Lorentz factor 7 . Z = 1 and 7 = 1 would correspond to a static 
pair-free medium. 


Z and 7 are functions of only one parameter of the front, which we now dehne. 
Let dE^/dS be the energy column density of the 7 -ray front [erg/cm^]. The GRB is likely 
beamed, yet it is convenient to dehne its isotropic equivalent = dnR^^dE^/dS). When 
an ambient electron is overtaken by the front, it scatters energy"^ Csc = aT(di?^/dS') = 
Era'llA tiR^. The relevant dimensionless parameter is 




^SC 

rrieC^ 




(3) 


Z{^) and 7 (.^) were found numerically in B02. They are calculated by solving the radiative 
transfer problem coupled to the dynamic problem of the medium acceleration. We briehy 
summarize the calculations here. 


Two processes play important roles in the 7 -ray transfer: Compton scattering 7 -|- 
7sc + and photon-photon absorption 7gc + 7 —> e"*" -t- e~ (here symbol 7 stands for a 
primary photon and 7 sc for a scattered photon). The same processes determine the deposited 
momentum and acceleration of the medium. The medium remains optically thin, so only 
single scattering is of interest, and the process 7 sc + e^ —> 7 sc + e^ is negligible. Nevertheless, 
there is a non-linearity in the problem because the scattering opacity of the medium is 
affected by creation and changes enormously across the front. 

The problem would be simpler if the scattered photons instantaneously converted to 
(Thompson & Madau 2000; Meszaros, Ramirez-Ruiz, & Rees 2001) — then there would 
be no need to solve the radiative transfer. However, the bulk of scattered photons never 
convert to and escape. The formal absorption free-path of a scattered photon in the 
radiation front turns out to be 3> A, which means that a scattered photon has a small 


^We assume that the shock wave is approximately adiabatic. Then Fghock ~ V^F, where F is the Lorentz 
factor of the postshock material (which we hereafter call “blast” for brevity), and I « R{1 — dshock) = 
^/2R\ock = ^/4rh 

"‘Accurate calculation of scattering includes the Klein-Nishina correction to Thomson cross section ut, 
however, in the definition of it is convenient to use ctt which is independent of photon energy. 
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Fig. 1.— Early stage of a GRB blast wave {R < i?ioad)- The short-burst (or thin-shell) regime 
is assumed: A < R/AT"^ (cf. eq. 1). The forward shock propagates in the e^-loaded and 
preaccelerated medium left behind the 7 -ray front. The pair-loading factor Z and Lorentz 
factor 7 of the preshock medium depend on the current radius R of the explosion (see 
Fig. 2). The inset schematically shows what happens in the 7 -ray front: some of the 7 -rays 
are scattered, lose their collimation and get absorbed by the collimated 7 -rays, producing 
e^. The created pairs are immediately Compton cooled in the 7 -ray front (B02), and the 
preshock medium is relatively cold {kT < rrieC^). The pairs are heated/accelerated when 
the forward shock reaches them, and the postshock produce the broad-band synchrotron 
emission. The whole structure — radiation front, e^-loaded preshock medium, blast, and 
ejecta — has a small thickness ~ R/T"^. 
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chance A/A^^ to create an pair. The exponential runaway of pair loading occurs when the 
smallness of this chance is compensated by a large number of scattered photons per electron, 
which requires a small scattering free path of the electron in the radiation front, A A. 
The length of exponential loading is then a = (AA^.^)^/^ < A. The characteristic loading 
radius i?ioad is dehned bya = A; A<^A<^ at this radius. 

At small R, the e^-loaded medium is accelerated to a signihcant Lorentz factor 7 (not 
to be confused with the photon symbol 7 in the reaction formulae). The photons scattered 
by the accelerating medium have the collimation angle 66 ~ 7 “^ and their chances to convert 
to pairs become completely negligible (a smaller angle between the scattered and primary 
photons implies a higher threshold for reaction 7 sc + 7 —^ e’*' + e~ and a smaller absorption 
opacity seen by the scattered photons). Therefore, loading at small radii is made by 
photons scattered in a small leading portion of the front (where the medium has not yet 
acquired 7 3> 1 ) and propagated across the front. 

The numerical solution to the transfer problem describes what exactly happens with 
the medium in the radiation front. The shell of ambient medium that is inside the front at 
a given moment of time has a certain velocity prohle P{w) where 0 < ci7 < A is the distance 
from the leading edge of the front (“entrance”). The medium velocity (3 increases from zero 
(7 = 1) at the entrance ci 7 = 0 to its hnal value at the exit zu = A. Lepton number per 
proton, Z, also increases from Z = 1 to its value behind the front. The front structure is 
described by the same functions 7 (.^) and Z(^) that describe the front evolution with radius 
if one substitutes ^ = zu/\ (B02). So, the front structure is described by the same functions 
at different times, i.e., it is self-similar. 

Since there is a gradient of the medium velocity inside the front — layers at larger zu 
move faster — one might think that caustics, i.e., internal shocks can develop in the front. It 
never happens at radii of interest. The radiation overtaking the medium dictates its velocity, 
and the medium has no time to develop a caustic because it quickly exits the front^ and is 
left behind with a uniform Lorentz factor 7 (the exit 7 gradually evolves on timescale R/c 
because ^ = A/X evolves as R~^ as the front expands). The medium left behind the radiation 
front is immediately picked up by the blast wave, on a timescale (Rlc){Thf « fl/c. 

The exact numerical solution for 7 (,^) and Z(^) is well approximated by a simplihed 
analytical model derived in B02 and summarized in the next section. 


^The exit time A/cy^ is much shorter than R/c at R > i?gap (§ 2.1). The fact that no caustics develop in 
the front at i? > Rgap formally follows from equation (12) in B02, which shows that density remains finite 
at all w, while a caustic would correspond to infinite density. At an earlier stage, R < Rgap, which is not 
considered in this paper, caustics do develop, see B02. 



2.1. Analytical description of e^-loading 


The exact Z(^) and 7 (.^) are approximated by the following analytical formulae (see 
eqs. 49, 62, 63 in B02; Z = n*/no in B02 notation). 


l[exp(^/6oad)+exp(-^/^ioad)] i ^ 'Caco 

= { {i/UcfZ,,, Uc <i< S^acc, 

3('C/'Cacc)-^acc ^ ^ Si^aco 

^ ^ ^ 'Caco 

7(0 = { (e/Occ)' Occ < e < S^acc, 

303(e/0cc)'/' e > S^acc, 


where ^acc = (5 + ln/ie) 6 oad, 


Zacc 


exp 


'Cload 


+ exp 


'Cload 


= 74/ie, 


( 4 ) 

( 5 ) 


( 6 ) 


and fie is the electron mean molecular weight of the ambient medium {fie = 1 for hydrogen 
and fie = for helium or heavy ions). The numerical value of .^load = 20 — 30 depends on the 
precise spectrum of the 7 -rays, however all the other relations remain the same. 

The standard GRB spectrum is a broken power-law with a peak at hu ~ nieC^, 


F. = Fi 


{hn/nieC^) hn < nieC^, 
(/iz//meC ^)““2 hn > nieC^- 


( 7 ) 


When observed from a redshift z ^ 1, the spectrum peaks at hu zz 250 keV, as reported by 
BATSE observations (Preece et ah 2000). .^load for such a spectrum was derived in B02, 


'Cload(cH) CI2) 


( 0:2 — Oil) 


(1 - q;i)(q;2 - 1) V20O 


0^2 — 0^1 


^KN 


1/2 


( 8 ) 


where ckn ~ 0.4 and 0=^2 -|- q; 2 ) In the examples below we £x ai = 0, so that 


A typical GRB has 0:2 


'Cload(o! 2 ) 


g\ 1/2 5^272^3/2 

7 ) a 2 -l 


(l + a 7 ''«. 


^1.5, which gives .^load ~ 120 . 


( 9 ) 


At sufficiently large radii, when the 7 -ray front has ^ < ^load, hs effect on the medium 
is negligible: the medium remains almost static (7 ~ 1) and e^-free {Z ^ 1). When the 
front has ^ > .^load, the runaway loading occurs. The number of loaded pairs depends 
exponentially on ^ as long as .^ < ^acc- At .^ > .^acc, the front acts as a relativistic accelerator 



9 


and the dependence of 7 and Z on ^ can be approximated by power laws (eqs. 4,5). The 
slopes of 7 (.^) and Z(^) change at ~ 3.^acc where Z ^ 10^ and the mass of injected is 
comparable to that of the ambient ions. An interesting effect takes place at .^ > .^gap ~ 3x 10^: 
then 7 (,^) exceeds the Lorentz factor of the ejecta Tej = 10^ — 10^. The radiation front with 
snch a high ^ pnshes the external medinm away from the fireball and opens a gap, so that 
the hreball moves in a cavity cleared by its own radiation front. 

The 7 -ray front expands with time, and its .^-parameter decreases as It starts 

at very high ( and then passes throngh ^ap, Cacc, and ^load at radii i?gap, i?acc, and i?ioad, 
respectively, 

cm, i?ioad = (5 -h acc5 -^gap -^acc /3. (10) 

Z(R) and j(R) of the preshock medinm are shown in Figure 2 for = 10^^ erg. They 
do not depend on the density of the ambient medium as long as the medium is optically thin 
and are entirely determined by the parameters of the radiation front. The figure shows the 
exact Z and 7 (cf. Fig. 1 and 3 in B02) and their analytical approximations (eqs. 4 and 5). 

Variations in GRB spectra have some effect on Z{R) and 'y{R) because they affect .^load- 
For instance, if 0:2 changes from 1.5 to 2 then .^load changes from 24 to 33 (eq. 9); the resulting 
changes in Z{R) and 'j{R) are shown in Figure 2. 


^acc 


7 X 10^® 


( E, 

\ lO^^erg 


2.2. Blast wave 

The prompt GRB emission can be produced quite early inside the hreball, preceding the 
development of the external blast wave at R ~ 10^^ — 10^^ cm. The millisecond variability 
observed in GRBs is often cited as an evidence for small radii of 7 -ray emission, R^ ~ 
10^^ — 10^^ cm. Regardless the precise R^, the 7 -rays are dehnitely emitted before the 
afterglow and then propagate ahead of the blast wave through the ambient medium. 

The three characteristic radii Rgap, Race and Rioad dehne four stages of the GRB explo¬ 
sion: 


I. R < Rgap. The hreball moves in a cavity cleared by its radiation front. The e^-rich 
external medium {Z > 10^) surfs ahead with 7 > Tej. 

II. Rgap < R < Race- The external blast wave has formed and the hreball sweeps up 

the e^-rich medium {Z ^ 10^) that moves with 1 7 < Fej. 
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R [cm] 


Fig. 2.— Pair-loading factor Z and Lorentz factor 7 acquired by a medium at a radius R when it 
is overtaken by a 7 -ray front. The front has isotropic energy = 10^^ erg. The exact numerical 
results are shown by solid curves and the analytical approximation (eqs. 4 and 5) by dashed curves. 
The characteristic radii Rgap, Race and i?ioad are indicated by vertical dotted lines. The condition 
7 = Tej defines Rgap] in this example Tej = 200 (shown by the horizontal dotted line). Z{R) and 
7 ( 1 ?) for different are found by simple re-scaling R {E^/10^^)^^‘^R. The GRB spectrum 
assumed in the calculation has the typical high-energy slope 02 = 1-5. To illustrate the dependence 
on 02 we also show Z{R) and 7 (i?) for a 2 = 2 (thin solid curves). 
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III. Race < R < -Rioad- The fireball sweeps up the static medium (7 ~ 1) dominated by 
e± (Z > 1). 

IV. R > -Rioad- The hreball sweeps up the static pair-free medium. 

The blast wave develops at R > Rgap, with its forward and reverse shocks. The preshock 
external medium is e^-rich and moving relativistically, which affects the jnmp conditions and 
the state of the postshock plasma. At R = Rgap the blast wave gently begins to sweep up 
the preaccelerated medium with a small relative Lorentz factor Trei ~ 1 (Tej ~ 7 ). With 
increasing R > Rgap, 7 falls off quickly and approaches unity at R = Race as 7 = (R/Racc)~®- 
Thus, the hreball suddenly “learns” that there is an interesting amonnt of slow material in 
its way and hits it with a large Trei. This resembles a collision with a wall and results in a 
sharp rise of the afterglow at R ~ Race- 

The explosion pictnre is modihed if the 7 -ray front is created by the shock wave itself 
(external model of the prompt GRB) rather than by internal dissipation in the hreball. Then 
a self-consistent blast-wave model exists without the gap. A sufficiently dense medinm can 
complicate that pictnre if it leads to an “electro-magnetic catastrophe” in the forward shock 
(Stern 2003) — a runaway prodnetion of 7 -rays, which may result in transient gap openings. 

The uncertainties in the mechanism of the prompt GRB play little role for the afterglow 
model developed below as we focus on the stage when the process of 7 -ray prodnetion is over. 
We consider R > Rgap, and the blast wave at these radii has practically no memory of the 
gap opening. We need only to know the parameters of the emitted 7 -ray front (eq. 7) which 
are taken from observations. Therefore, the following calculations apply to both internal 
and external models of the prompt GRB emission. 


3. Afterglow calculation 

The production of afterglows without the e^-loading effects, i.e. with Z = 1 and 7 = 1, 
was stndied previously in detail. The electrons were assumed to receive a signiheant fraction 
Ce ~ 0.1 of the shock energy (the remaining fraction 1 — Cg is carried by the postshock ions). 
This leads to a high energy per electron, and the early high-T blast wave emits hard X-rays 
in the fast-cooling regime: the electrons immediately radiate the received energy, practically 
at the same radius where they are shocked. 

The loading changes the picture: a lot of leptons now share the energy received from 
the shock. As a resnlt, the mean energy per postshock is reduced by orders of magnitude 
and emit much softer radiation. Furthermore, the postshock can be in the slow-cooling 
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regime so that their emission remains signihcant at i? > -Rioad where no new pairs are added 
to the blast wave. 


3.1. Formulation of the problem 

We focus on a sufficiently early stage when F > 9~^ where 9 is the opening angle of the 
explosion; then the expanding jet behaves like a portion of a spherically-symmetric explosion 
and we can neglect the beaming effects. Three further technical assumptions simplify our 
calculations: 

1. — The shocked ambient material is assumed to have a common Lorentz factor F 
(relative motions in the postshock material are subsonic and the assumption of common F 
is good within a factor of l/v^). 

2. - The postshock material is assumed to be in pressure equilibrium, i.e. P = const 
throughout this material. 

3. — When calculating the observed luminosity we neglect the finite thickness of the 
blast wave (= the distance between the contact discontinuity and the forward shock). The 
blast-wave radiation is then characterized by its total instantaneous luminosity emitted 
at a current radius. Inclusion of time delays between photons from different subshells of the 
blast material would change the result by a factor of a few. A precise calculation of this 
factor would require to relax assumptions 1 and 2. 

These assumptions will allow us to derive the instantaneous luminosity L^, at a given 
frequency z/obs- We are interested in the soft emission here, from IR to soft X-rays, and 
especially optical, z/obs ~ 0.5 x 10^^ Hz. 

Consider the moment of time when the blast wave reaches a given radius R. The total 
postshock ambient mass at this moment is® 

m{R) = / 4nR‘^po{R)dR. (11) 

Jo 

Each subshell Sm of this material has its own history. It used to be a shell of ambient 
medium located at a radius R < R, which was loaded with pairs and preaccelerated by 
the prompt y-rays, then immediately shocked and picked up by the relativistic blast wave. 

^Hereafter quantities taken at the current radius R are marked by tilde to distinguish them from the 
corresponding quantities at radius R(m) where shell (m, m -I- 5m) was shocked. 
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The radius where all that happened is related to m by equation 

l‘R{m) 

m= 47r-R^Po(-R)d-R. (12) 

Jo 

Then Sm cooled radiatively and adiabatically as the blast wave expanded to the present 
radius R. 

We need to evaluate the contribution of each Sm to the current luminosity L^^R) and 
therefore have to resolve the blast wave material in its (Lagrangian) mass coordinate 0 < 
m < fh, even though we do not resolve it in radius and assume that all Sm are located at 
the same R. m —0 is the material which was swept up hrst (at small radii), and m ^ m 
is the currently shocked material. L^, is the sum of current emissions SL^ from all Sm. 

In afterglow models without loading the emission from m m is negligible and the 
instantaneous luminosity is dominated by the recently shocked mass shell m ^ fh. A special 
feature of pair-loaded blast waves is that their peaks at m m. Therefore an accurate 
integration over m is needed. 

To calculate SLi, from each Sm we need to know the current spectrum /( 7 e) and 
magnetic held B in Sm. 


3.2. Magnetic field 

The magnetic held in shell Sm in the blast wave of radius R is given in terms of the 
customary equipartition parameter eBim) < 1 , 

UB{m) = = eB{m)U, (13) 

Stt 

where U = 3P is the current energy density in the blast wave and P is its pressure. P and 
U are assumed to be constant throughout the postshock material in our simplihed model. 
They are functions of R only and related to the ambient density by the jump condition at 
the shock front, 

U = 3P = dpoc'freif, (14) 

where 

fj-ei = ^ ^ ~ (when 7 < f) (15) 

7(1 + W 

is the current Lorentz factor of the preshock medium relative to the blast wave, and 7 = 
(1 — ig the Lorentz factor of the preshock medium in the lab frame. The jump 
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condition (14) takes into account the compression of the ambient medium preaccelerated to 
a Lorentz factor 7 (Madau & Thompson 2000; B02). Thus, the magnetic held can be written 
as 


B = f 

3277 Cb PqC^ 

1/2 

= 0.39 f 

eghoPe 


- 7(1+ /3) . 



. 7(1 + /3). 


where ho is expressed in cm 


G, 


(16) 


€ 5 ( 771 ) can change with R\ the magnetic held evolves as 5m expands from the radius R 
where it was shocked to the current radius R. The toroidal held component is dominant in 
the expanding shocked plasma and, if the magnetic hux is conserved, B evolves as 


/ ~ \ 3/4 ^ 
B{m) p R (^ 

B{m) p R \P j -R’ 


(17) 


where p is the proper mass density of baryons and we have used the adiabatic equation of 
state P cx Using equation (14), we can rewrite the hux-conservation condition as 


B{m) _ f fr^elPo A ^ R ^ _ f r^relPo ^ ^ f R\ 

B{m) yTTrelPoy -R’ es ^TTrelPoy y-Ry 


(18) 


It shows that €5 grows downstream of the shock. This growth, however, cannot proceed 
beyond equipartition: €5 then saturates near unity. 

If the hux is not conserved (the magnetic held may reconnect/dissipate in the blast) the 
evolution of postshock B is diherent. Therefore, we will not specify €5 in the calculations 
until the last step when we give examples. Then, besides equation (18) we will also consider 
the case of €5 = €5 = const, which is a second reasonable prescription for the magnetic held 
evolution. 


3.3. Distribution function of electrons/positrons 


The customary phenomenological shock model of GRBs assumes that the electrons (and 
positrons) are impulsively accelerated at the shock front with a power-law spectrum. 


/(7e) 


0) 7e ^ Imi 

{'Je/1 m) j 7e 7 7m. 


(19) 


This initial spectrum is injected in 5m at the radius R{m) where 5m is shocked; the Lorentz 
factor 7 e of the accelerated leptons is measured in the huid frame. The total injected energy 
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of nonthermal in this frame is 


6E= JemeC^fije) dje = eeTj-elSmC^, 


( 20 ) 


and the total number of nonthermal leptons in Sm is 


5N= /(7e) d7e = 


6m 


z. 


( 21 ) 


Here p > 2 and eg < 1 are phenomenological parameters of the electron acceleration. 
From equations (20) and (21) one hnds 7 ^ and K, 

hrel ee(p - 2)mp 


'Im 


(Z//ie) (p-l)me 

K = {p-l) — . 


( 22 ) 

(23) 


When a shocked 6m expands from Rio R the injected nonthermal spectrum is modihed 
by two effects: 


1. — Adiabatic cooling shifts the whole distribution as 7 e oc retaining the 

power-law shape. The minimum Lorentz factor of the nonthermal spectrum, 7 ^, changes by 
the factor 


A 


p 

1/4 

Poh relT 

p 


Poh relT 


(24) 


7,n(m, R) = R). 


(25) 


and the normalization K of the spectrum changes as K = K/A. 


2. — Radiative cooling cuts off the distribution at high 7 e. In most of the models 
considered below, the exact position of the cutoff is not important because we focus on the 
low-frequency radiation. The cutoff can be estimated as follows. 

The radiative cooling of 6m can peak at any radius R' between R and P, and then 
the cutoff is shaped at this radius. Electrons with Lorentz factor 7 e are cooled with a rate 
(assuming isotropic pitch-angle distribution). 


7e 


4 (Tx 
3 rrieC 


r )/2 


7e^ 


(26) 


where B' = B{m,R') and Pj = Us{m,R') is the energy density of soft (synchrotron) radi¬ 
ation in 6m at radius R'. We assume that the bulk of synchrotron radiation in the fluid 
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frame satisfies /ir'fiuid < ^ef^he and scatters off with Thomson cross section. Using 
eqnations (13,14) we get 

7. = -15—<TTc4r;.,r'(i + cy,l ( 27 ) 

o TTh^ 

Here C = U'g/U'^ is the relative contribntion of inverse Compton scattering to the cooling 
rate. 


The characteristic cooling Lorentz factor 7 ' at R' is dehned by the condition that the 
cooling timescale 7 e/| 7 e| eqnals the expansion timescale {R! — R)/V'c. This yields 


7'(m) 


Srrie 


(28) 


Note that C depends on the radiated energy U', which in tnrn depends on 7 '; therefore 
eqnation (28) is implicit. It can be solved as follows. The total radiation density in 5m 
is comparable to the lost energy: ~ {'l'ch'rn)'^~^^eU' if 7 c > 7 ^ and ~ eJJ' 

otherwise. Its synchrotron fraction is U'g/U'^^ = U'^/{Ug + U'^). This allows one to express 
C = Ug/Ug in terms of 7 ', 


C\C' + 1) 


^'b 


1 7 ^ < 7 ^ 

ic ^ im 

ic > im 


(29) 


Now we have two eqnations (28) and (29) which can be solved for C and 7 '. 


This estimate assnmes that Ug{m) is prodnced locally at a given m and neglects the 
transport of synchrotron radiation across the blast wave. Inclusion of transport further com¬ 
plicates the calculation of C'{m) and 7 c(m). In this paper, we avoid models with signihcant 
Compton cooling and use the simple estimate of 7 ' (eq. 28) with C = 0. The consistency of 
this estimate requires ^ 1, which can be checked using equation (29). An estimate of 

Compton cooling that includes the radiation transport will be given in § 6.3. 

The distribution cutoff % in 5m at the current radius R is 


7c = 7^', A' = 


PoTrelT 

PoKel^' 


1/4 


(30) 


It is shaped at the radius R' where ( 7 ' A') is minimum with 7 ' calculated as described above, 
and its evolution from R' to R is determined by the adiabatic cooling factor A'. 

If 7 m < 7 c (slow-cooling regime), the current nonthermal distribution in 5m is a 
power-law with slope p and normalization K, extending from 7 ^ to 7 c. If 7 c < 7 m (fast- 
cooling regime), all nonthermal leptons 5N pile up near a single Lorentz factor 7 e ^ 7 c. In 
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any case, the number of nonthermal in Sm is given by equation ( 21 ), so their distribution 
function at R can be written as 


/(7e) =SN X 


(7e/7m) ^ - 7m) hf(7c “ 7e), 7m < 7c- 

*^(7e 7c)) 7m 7 7c, 


(31) 


where is Heaviside step function and 5(...) is Dirac ^-function. 


3.4. Self-absorption 


Self-absorption can affect synchrotron emission from with low 7 e. The low- 7 e popu¬ 
lation is created at the very beginning of pair-loading {R < -Race) when the preacceleration 
7 is signiheant, the forward shock is relatively mild (Trei T), and the pair-loading factor 
Z is comparable to 10^ (cf. eq. 22). The low-energy postshock are slowly cooling and 
produce low-frequency radiation that is self-absorbed. 


Self-absorption is, however, not important for blast-wave radiation with frequencies 
u ^ 10^^ Hz that we consider thereafter. This radiation is dominated by the blast-wave 
material with sufficiently high 7 ^ so that self-absorption can be neglected for all 7 e > 7 m- 
The condition of small self-absorption can be written as Um < 7 m^^eC^(z^/c)^, where Um is 
the density of synchrotron radiation produced by with 7 e = '-fm, and we assume thereafter 


R 

Um ~ ImUsCTT — < 'JmmeC^ 


uy 


(32) 


This roughly corresponds to z/ > 10^"^ Hz. 


3.5. Synchrotron luminosity 

Given the current magnetic held B and nonthermal spectrum /( 7 e) in Sm, it is 
straightforward to evaluate its contribution SLi, to the instantaneous luminosity of the blast 
wave. 

The synchrotron spectrum of with Lorentz factors 7 e in the huid frame peaks at the 
frequency (assuming isotropic pitch-angle distribution) 

Z/flmd^0.15^7e- (33) 

rrieC 

The intensity-weighted Doppler shift to the lab frame is given by 

— ^ r 
^ 2 -L 


( 34 ) 



18 


and the corresponding observed freqnency is 


V 


^obs — 


1 + ;^’ 


(35) 


where is the cosmological redshift of the bnrst. From eqnation (33) we hnd 7 e of whose 
synchrotron spectrnm peaks at a given v. We denote this characteristic Lorentz factor by 
W, 

5 Z/TTlpCl 


Iv = ly{m,R) 


1 1/2 


TeB 


(36) 


The emitted synchrotron power by the electron is Eg = and the syn¬ 

chrotron Inminosity from an popnlation with distribntion dA^/dye is approximately given 
by 


dL 


dL 


^ d log V 2d log 7 i, 


- , diV\ . , , 


(37) 


7e=7i- 


The popnlation in bm peaks near if 7c > 7m (slow-cooling regime) and near 7 c if 
7 c < 7m (fast-cooling regime). In the first case, the synchrotron Inminosity of 6m can be 
written as 




X 


W < 7m < 7c, 
(w/7m)^"^ lm<% < 7c, 
0 7m 7 7c 7;y, 


where 


5U 


aTcB'^^^ 5 meC^a^B 

- 7 ,, 0 A/ =-^ V oN. 

127r 1271 e F 


(38) 


(39) 


Note that 5{vLy) is originally dehned as energy emitted per nnit time in the hxed lab frame, 
however, we calcnlate it in the flnid frame nsing the Lorentz invariancy of vL^. is a 

maximnm Inminosity that would be emitted at v at the most favorable condition y^ = % [in 
this case 25{i'L^) equals the energy loss rate of the dominant population with yg ~ 7m]- 
Similarly, in the fast-cooling regime, we have 


6{uL^) = X 


(Ru! 7c) ^ 7^- 7 7c 7m, 

0 7i. > 7c. 


(40) 


The nonthermal with Lorentz factors 7 e = 7m and 7 e = 7c emit radiation at frequencies 


~ eB 

z/m = 0.2F-7^, 

mpC 


z>c = 0 . 2 F-^ 72 , 
mpC 


( 41 ) 
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and the spectral luminosity 5Ly/5m at a Lagrangian coordinate m can be written as 


6L, 

6m 


( 


6L^^ 

5m 


X < 




0 


V <Vra< Vc, 
Vm<y < Z>c, 

V <Vc< Vm, 

V > Z>c. 


(42) 


An explicit formula for 6U^^^ = 6U^^^/v is found using equation (21) for 6N and equa¬ 
tion (16) for 5, 


^^max 

5 

meC^ (Tt B 5N 

6m 


e V 

5m 


30 

IbUo 

-he7(l + /3). 

1/2 


327r' 


rrip 


esuorl 


f^emp 7(1 + p) 


1/2 


(43) 


where ^ jm^f? ~ 2.81 x 10 ^^cm is the classical electron radius. Note that 6Up^^j6m 
does not depend on v. 

The total instantaneous luminosity of the blast wave is found by integrating over m, 


U{R)= r^-^Sm. 

Jo Sm 


(44) 


3.6. Observed flux 


Observer at a distance D much smaller than the Hubble scale would measure the spectral 
luminosity (e.g. Rybicki & Lightman 1979), 

L,, (45) 

and the spectral flux For cosmologically large D, the flux formula is 

modified by two effects: D is replaced by the luminosity distance and the observed frequency 
of radiation is redshifted, z/obs = (1 + 


where 


_ f\l + z) 
47rz/obs-D^ 37iD‘^ 


Fq 


(i + z- VTTz^ ^ 2.6 X 10^® (i + z- VTTz^ 


cm, 


(46) 

(47) 
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and Hq k, 70 km/s/Mpc is Hubble constant. The observed optical magnitude in V or R 
band is related to by 

my = 8.873 -2.5 IgF^, m^ = 8.645 - 2.5 IgF^, (48) 

where T7 is in units of Jy = 10“^^ erg cm“^. 

The flux T7(i?) is received at the observer time 

tohs{R) - ^ {I + z) (49) 

21 

after the arrival of hrst photons (prompt y-rays) from the explosion. More exactly, it is 
received during a time interval Afobs ~ (i?/2r^c)(l + z) because of the spherical curvature 
and hnite thickness of the blast wave, and equation (49) gives only an approximate arrival 
time. The numerical factor in this equation can vary around 1/2 by a factor ~ 2, depending 
on the blast wave dynamics r(R). We hereafter use the approximate equation (49) with the 
hxed factor 1/2. This degree of accuracy is consistent with our approximate treatment of 
the blast wave as a constant-pressure shell. 


4. Numerical examples 

It is straightforward to calculate the afterglow emission F^(tohs) using the formulae of 
§ 3 and taking the integral (44) numerically. A simple illustrative afterglow model has the 
following parameters: 

1. Isotropic energy = 10^^ erg and initial Lorentz factor Tq = 200 (the reverse shock 
is assumed to be non-relativistic and Tq ~ Tej). 

2. External density no(R) = const and He = 1 (uniform hydrogen medium). We will 
consider no = 0.1, 1, 10, 10^, and 10^ cm“^. 

3. Parameters of the postshock leptons: Ce = 0.1 and p = 2.5. Most of the blast-wave 
energy is carried by relativistically hot ions which do not exchange energy with the electrons 
and the blast wave is approximately adiabatic. 

4. Magnetic equipartition parameter cb immediately behind the shock. It will be taken 
equal to 10“®, 10“^ or 10“^ in the example models. eB{ni,R) can evolve in the postshock 
material of the blast wave as it expands and decelerates. We will consider two cases: magnetic 
flux is conserved in the postshock region (eq. 18) or 6 ^( 771 , i?) = cb is constant in all shells 
m at all times. 
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5. The 7 -ray front is described by its isotropic energy In the examples below we 
assume = i?ej = 10 ^^ erg, i.e., half of the total explosion energy is initially emitted in 
the prompt GRB. We also assume the standard GRB spectrum with a high-energy slope 
0:2 = 1.5 which gives ^acc = 120 (§ 2). 

The blast wave has a constant T ^ Tq until it approaches the characteristic deceleration 
radius i?dec where the swept-up mass m satisfies Tgmc^ = E^y At larger radii r(i?) quickly 
approaches the self-similar solution of Blandford & Mckee (1976): T = (17/12)^/^ro(R/Rdec)~^^^. 
We will use a simple approximation: r(i?) = Tq at R < Rdec and r(R) = ro(R/Rdec)~^'^^ at 
R > Rdec- It has sufficient accuracy, adequate to our simplihed hydrodynamical model (§ 3). 

The results of numerical calculations are shown in Figures 3-7. In all examples, a 
cosmological redshift 2 : = 1 is assumed. 

Figure 3 shows the optical light curve Eyitohs) found at no = 10 cm“^ and = 10“"^. 

It has two peaks. The first peak is a result of loading (for comparison we show the 
results with Z = 1 and 7 = 1, i.e., with neglected impact of the 7 -ray front on the external 
medium). The figure compares the results obtained with magnetic flux conservation and 
eB{m,R) = const. The emission is stronger if the magnetic flux is conserved because it 
implies a higher cb] then the peak of radiation is reached at R = Rdec (see § 5). 

The second peak is produced by the pair-free part of the blast. It is well described by 
the standard afterglow model and corresponds to = n (see § 5.1). The time of the second 
peak fpeak oc 

Figure 4 shows the corresponding light curves in the infrared and X-ray bands. The 
relative contribution of to the afterglow is higher at lower frequencies, and their effect 
is negligible in the hard X-ray band. Gomparing the light curves in IR, optical, and soft 
X-rays, one can see that the spectral flux of the peak is comparable in all three cases, 
indicating a small spectral index of radiation. 

Figure 5 shows the instantaneous broad-band spectrum emitted by the blast wave at the 
deceleration radius Rdec (which is 4.6 larger than R^cc in the example model). The small e^- 
dominated fraction of the swept-up mass dominates the emission in soft bands. Its spectral 
slope is small: |q;| ^ 0 . 2 . At high frequencies, z/obs > 10^^ Hz, the spectrum is dominated 
by the pair-free material behind the current position of the shock and well described by the 
usual pair-free model. In particular, the standard breaks at 10^® Hz u 

and z/c ~ 10^'^ Hz seen. [In the upper panel, the second (cooling) 

break occurs at a smaller frequency ~ 10^® Hz and is smooth because happens to be 
close to Pm-] The spectrum of the emission component has a break at z/ ^m(^load 7 -R)? 
which is ~ 10^^ Hz in the upper panel and ~ 3 x 10^^ Hz in the lower panel, and cuts off 
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Fig. 3.— Example of the optical light curve in the R-band, i^obs = 5.45 x 10^^ Hz, for a 
’’canonical” GRB explosion: adiabatic blast wave with Eej = = 10^^ erg, Fq = 200, 

z = 1, 77-0 = 10 cm^, eB = 10“^, = 0.1, p = 2.5. Left axis indicates the observed flux in 

mJy and right axis — the corresponding R-magnitude. Calculations shown by solid curves 
assume conservation of the postshock magnetic flux. Dashed curves show the results with 
eB{m,R) = const = cb assumption. Thin light curves would be obtained if the pairs were 
neglected (Z = 1 and 7 = 1 ). 
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tobs [s] 


Fig. 4.— Light curves at i^obs = 10^^ Hz (upper panel) and i^obs = 10^^ Hz (lower panel) for 
the same model as in Figure 3. The meaning of line types is the same as in Figure 3. 
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Fig. 5.— Instantaneous spectrum emitted by the blast wave at i? = i?dec is shown by the 
solid curve (the model parameters are the same as in Fig. 3). Dashed curve shows the 
contribution of the e^-loaded shell m < mioad- Dotted curve shows the spectrum that would 
be found in the absence of y-ray front [Z = 1 and 7 = 1); it is approximately equal to the 
spectrum of pair-free material m > mioad- Upper panel: conservation of postshock magnetic 
flux is assumed. Lower panel; eB{m,R) = const = cb is assumed. 
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Fig. 6. — Optical light curves found for external densities Uq = 0.1, 1, 10, 10^, and 10^ cm“^. 
All other parameters are the same as in Figure 3. Calculations shown by solid curves assume 
conservation of the postshock magnetic flux. Dotted curves show the results with esim, R) = 
const = cb- 
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Fig. 7.— Optical light curves found for = 10“®, 10“^, and 10“^. All other parameters 
are the same as in Figure 3. Calculations shown by solid curves assume conservation of the 
postshock magnetic flux. Dotted curves show the results with esim, R) = const = cb- 
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at u > Uc{m\oaAi R)■ The cutoff frequency is smaller in the upper panel because the flux 
conservation gives a stronger magnetic held. 

Figure 6 shows the optical light curves for different Uq and hxed = 10“^. It also 
illustrates the effect of magnetic hux conservation by comparing it with the results obtained 
at es(m,i?) = const. The emission component quickly increases with increasing no. It 
scales as as we show in § 5, and this scaling is observed in Figure 6 for constant-ee 
models. 

Figure 7 shows how the optical light curve depends on the shock parameter e^. At 
cb = 10“^ the are in the fast-cooling regime (z/c < r'm) and most of the pair energy is 
emitted at earlier times. This produces a bump in the light curve before or at the deceleration 
radius, depending on the behavior of magnetic held downstream of the shock. 

5. Analytical calculation 

In § 5.1 we derive the luminosity of a relativistic blast wave without a y-ray precursor, 
when the e^-loading and preacceleration ehects are neglected {Z = 1 and 7 = 1). The 
pair-free blast wave was studied in a number of previous works (e.g. Sari, Piran, & Narayan 
1998; Granot, Piran, & Sari 1999; Panaitescu & Kumar 2000) and used extensively to ht 
observed afterglows. We show that the previous results can be obtained with a diherent 
technique where the produced luminosity is calculated as an integral over the Lagrangian 
mass coordinate m (eq. 44). This approach is extented to e^-loaded blast waves in § 5.2. 

Where the ambient density prohle needs to be specihed we assume a power-law, so that 
the swept-up mass is 

R y 

-^acc / 

k = 3 describes a uniform medium no{R) = const and k = 1 describes a wind-type medium 
no oc R~^. The derivation in §§ 5.1 and 5.2 applies to both cases. In § 6 we consider in more 
detail the case k = 3. 





5.1. Pair-free afterglow 


From § 3.4 we have at Z = 1 and 7 = 1, 
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> dm. 


( 51 ) 


Symbol ”0” in marks the neglect of loading. Three qnantities in the integrand depend 
on the Lagrangian coordinate m: e^, 7m, and 7 c. Qnantities eBijn) and Vmijn) can vary 
signihcantly in the region m m, however, this region makes a negligible contribntion to 
the integral. Therefore, one can approximate and z>m as constants: eBijn) ~ esim) and 
Umi'fn) ~ ^m(hr). Only i>c(m) cannot be assnmed constant: i>c(m) oo at m ^ fh. 


A characteristic Uc can be dehned at m = m/2: it will represent the cooling cntoff in 
the main part of the postshock material. The value of 7c(m/2) is shaped as the blast wave 
expands from R = (the shock radius of m/2) to the current radius R. It is given by 

(see § 3.3), 


7c(m/2) 


a mg 

(T'j'Tnp /ig ei^eB^ R 


( 52 ) 


where a ~ (3/16)(l — 2 ^/^) ^ is a numerical factor. We assume here a small Compton factor 
C < 1 for simplicity. 


Let us dehne = jyimjZ) ^ Jvim), 7m = jmim/2) ^ 7 m(Ai), and 7 c = 7c(m/2); 
they are given by equations (36), (25), and (52), respectively. The numerical values of these 
Lorentz factors and the corresponding synchrotron frequencies are^ 


7. = 2.7x lO^z/Jf r2-'(esno/ig)-'/", (53) 

7 c = 8.2 X 10^ai?/g^r7^ (eBno/ig)“\ z/c = 9.2 x 10^® (esUo/ig)"^'^^ i?))/Hz, (54) 

7m = 1.84 X 10 ^r2i//ig, z/m = 4.6 X 10 ^°r2 (eBno)^'^We'^^ Hz, (55) 


where 


^ejp- 
iv-i) ' 


(56) 


There are two possible cases: 


^Hereafter in § 5.1 we drop tilde for R, is, and F because there is no need to consider m fh for the 
pair-free afterglow and distinguish between the postshock and current quantities. 
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Case I. h>c > Most of the blast material contributes to the observed luminosity at the 
frequency v. The luminosity integral over m is then estimated as k, m{5L/5m) where 
5L/5m is taken at m = m/2 and given by equation (51) with z/^ and Um replaced by and 

^m- 


Case II. z/c < z/. Most of the blast material does not emit any synchrotron radiation at z/ 
because its spectrum is cut off at a lower frequency Uc- Only a mass fraction (m — nic) m 
just behind the shock front contributes to L^, where nic is dehned by Vcij^c) = The 
condition z/^ > z/ is satished in the mass interval me < m < fh and the luminosity integral 
over m can be estimated as uL^ ^ {m — mc){5L/5m). 

The cooling Lorentz factor 7 c(m) increases toward the shock front m as 7 c(m) oc (i? — 
R)~^ oc and we have 


7c(m/2) 1 — / m 

7c(m) 1 — (1/2)^/^ V m 


(57) 


where so dehned b{m) is a slowly varying function: 5 = 2 at m = m/2 and b = [A;(l—0.5^/^)] ^ 
at m —m. Equating 7 c(m) and we hnd nic, 


m-nic _ ^ 
fh Iriy 


1 

b\v) ’ 


(58) 


with 6 = 2 at z/c = z/ and b = [k{l — 0.5^^^)] ^ at z/c z/. One can take 6 2 in all cases. 


The results can be summarized as follows, 

1/2 
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(59) 


(60) 


It agrees with the usual afterglow description (e.g. Sari et ah 1998). 

After the deceleration radius, we have for an adiabatic blast wave T^mc^ = Egj and then 


FAR) = 


Eej(l + ;^) 
2tiD‘^c^ 6m 
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One can see from this equation that the observed flux at a fixed v scales as ng' It can be 
found as a function of observer time fobs using the approximate relation fobs ~ (-R/ 2 r^c)(l+ 2 :) 
and the deceleration law for an adiabatic blast wave, which gives 


^obs 


Rmc 

2Ee,\ 


( 1 +^). 


(62) 


5.2. Pair-loaded afterglow 

We now calculate the instantaneous luminosity of the blast wave L^{R) taking into 
account loading. The pairs dominate the material with Lagrangian coordinate m < mioad- 
The luminosity from this material can be written as 

/■"ii 

L^{R) = / mi = min{mioad,br}. (63) 

Jo 

At radii R < i?ioad (hr < mioad) the whole blast is dominated and its total luminosity 
Li, = L^. At radii R > i?ioad ijh > mioad) the total luminosity is a sum of two parts: 
from small m < mioad and luminosity from the e^-free material mioad < m < fh. The 
latter peaks near m (§ 5.1). 

The total luminosity can be written as the sum, 

U = L^ + Ll. (64) 

At small radii R < R\oaA we have just one integral however, may be formally 

kept in equation (64) and interpreted as a luminosity that would be obtained at Z = 7 = 1. 
It would not change Ly because Lj)" at i? < i?ioad in the soft bands of interest. 


5.2.1. Calculation of 

Pair-loaded blast waves have a special feature: the peak synchrotron frequency r'm(nr, R) 
varies enormously with the Lagrangian coordinate m at m < mioad (see an example in 
Figure 8 ). The peak frequency is given by 

~ eB f^r fin 

hm(nr) = Umirn, R) = 0.2r- 7 ^ = 4.6 x 10^^-- - - y - yjy Hz, (65) 

^eC ^ 7(l + /3)[7(l + /3)]3/2Z2ny' 

where we have taken into account the adiabatic cooling factor A (§ 3.2) and assumed a 
slow radiative cooling of e^. The factor appearing in this expression is a very 
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£ 0.1 1 10 100 0.01 0.1 1 10 100 

/ acc / acc 

Fig. 8 .— Instantaneous blast-wave emission as a function of Lagrangian coordinate m at radius 
R = IZdec- The blast-wave model is the same as in Figure 3. The results that would be obtained 
without the 7 -ray front effect (Z = 1 and 7 = 1 ) are shown by dashed curves. Upper panel: 
Observed synchrotron peak frequency r'm(m,i?dec) and cutoff frequency r'c(m, ildec) corrected for 
cosmological redshift (i.e. divided by 1 + z). The horizontal dotted line indicates the i?-band 
frequency t'obs = 5.45 x 10^^ Hz. Lower panel: Distribution of the produced optical flux over 
m. The emission from the e^-loaded mass m < mioad sharply peaks near mace, at a specific 
m* where + z) = t'obs- The pair-free emission component peaks at the current position 

of the shock m = m = rridec- Left: Conservation of postshock magnetic flux is assumed. Right: 
€B{m,R) = const = cb is assumed. 
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steep function of m, R{m), or ^{R), whichever is chosen as a Lagrangian coordinate in the 
blast wave. It varies by five orders of magnitude when ^ varies from .^acc/2 to 2,^acc; which 
corresponds to a narrow range of R from i?acc/v^ to v^i?acc- 


The steep variation of with m implies that a narrow mass shell Am ~ m^cc dominates 
in a broad range of u. This fact enables a simple estimate of To a first approximation, 
the number of whose emission peaks at a given z/ does not depend on z/ and equals the 
number of particles in the shell Am ~ mace, 

iV*~^acc^ = 74^. (66) 

fJjeTTlp Tflp 


The produced luminosity dL^/dlogz/ approximately equals EN^/2, which gives the spectral 
luminosity, 


L^iR) = 


5 mgC^UT B: 
127r e 


iV, = 30 f 
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•^acc^acc fut 
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(67) 


Here Hacc is the current magnetic field in the shell Am and e^acc is the corresponding 
magnetic parameter divided by the current energy density of the blast wave). In 

the first approximation, = 1, i.e. does not depend on u and the spectral index of 
emission a ~ 0. 


Calculation of by accurate integration over m gives the correction factor f^, and 
confirms that ~ 1 in a broad range of frequencies. The correction factor is derived 
analytically in Appendix and shown in Figure 9. It is conveniently expressed as a function 
of z//z>acc where z>acc = t'acc(7?) = VmiRiacc) R) (thus, z>acc describes the peak of synchrotron 
emission from material shocked at i? = i?acc after it expanded to the current R), 
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( 68 ) 


This frequency is near the optical band. 


The basic reason for the small slope a is the steep dependence of the pair-loading process 
on radius (the big power 17/2 appears in the problem). 


5.2.2. contribution to the observed luminosity 
The total observed flux of the blast-wave radiation is (see eqs. 46 and 64) 


= F/ + F/ = (LJ + Ll) , 


(69) 
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Fig. 9.— Correction factor fiy (with is* = eeacc, see Appendix). Solid curve shows the 
case of k = 3 (uniform medium), and dashed curve — k = 1 (wind-type medium). The 
characteristic frequency z>acc is given by equation (68). 
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where the pair-dominated part is given by equation (67) and the pair-free part by equa¬ 
tion (59). This analytical result is compared with the numerical calculations in Figure 10. 
The analytical calculation of is more accurate because the integral peaks in a narrow and 
well dehned region of m. The standard pair-free part is still reasonably well approximated 
and we can compare the two components analytically. 

Using equations (59) and (67), we get 

= 74 ^. ( 70 ) 

m Qu \ eB J 

peaks in the X-ray band during the early afterglow, and 1 in IR, optical, and UV 

bands. By contrast, has the spectral factor ~ 1 in a broad range of soft frequencies 
from IR to soft X-rays (§ 5.2.2). 

The e^-loaded material dominates the observed flux in the soft bands for two reasons: 
(1) the number of emitting particles is increased by the factor Z ~ Z^^ = 74)Ue and (2) their 
synchrotron peak frequency is shifted to soft bands (because the mean energy per postshock 
particle is reduced). In addition, e^acc > if the postshock magnetic flux is conserved in 
the blast wave (see § 3.2). As a result, F^ 3> even when the blast wave expands to 
R > Rioad where the e^-loaded fraction of the blast material is very small, mioad ^ hr. 


6. Blast wave in a uniform medium 


We now consider in detail the case of a uniform hydrogen medium: no(i?) = const 
{k = 3) with )Ue = 1, and aim to hnd the observed flux F^ = F^ -|- as a function of fobs 
and z/obs- We shall assume i?ioad < Rdec where 
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(71) 


We shall also assume that the blast wave is approximately adiabatic and its Lorentz factor 
is given by 


f = rox 


1 R ^ R(^qc} 

(R/Rdec)"^/^ R > Rdec- 


(72) 


The approximate relation between the observer time fobs, the blast wave radius R, and its 
Lorentz factor T reads (see § 4), 


R(fobs) ~ 2f^c 


fobs 


1 + ;^ 


^^0*^(^obs/l T z') fobs ^ ^dec, 
Rdec(Rbs/^dec) ^ Rbs ^ ^dec, 


( 73 ) 
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Fig. 10.— Optical light curve calculated numerically in Figure 3 is here compared with the 
two components and calculated analytically (dashed lines). 
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where 


tdec = (1 + ^)^ ~ 194 (1 + ;^) 
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is the deceleration time. We will use below the following explicit relations at i? > i?dec, 
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(76) 


If the postshock magnetic flux is conserved in the expanding blast wave, the evolution 
of 63 is described by equation (18), 
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(77) 


where is the magnetic parameter immediately behind the shock when the blast wave had 
radius R{m), and j3{m) is the preshock medium velocity at that moment. 


6.1. Pair-free part 


The e^-free afterglow is described by equations (46) and (59), 
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We focus here on low frequencies u < Uc, Vm and keep only the first two lines in the general 
expression for (eq. 60). Substituting equations (55) and (54) for Um and Vc and using 
a ~ 1 (§ 5.1), one gets 


= 0.013 (e^no) z/< < z/^, 

= 0.48 z/^f .Rif (eijno)^/^ v < 


(80) 


Note that z/ = (1 + z)i'ohs in all the formulae. Substitution of /?(tobs) and f (fobs) gives as 
a function of fobs- 
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In particular, at tohs < tdec we have F = Fq = const and get 

f y/»/ io-» ,4 /s4/3,j-2/3(i 

Jy Dl, \ 10 -* nl,s (1 + 


u < Um < nc, 
z/ < z/c < z/„. 


(81) 


At tohs > tdec, we have for an adiabatic blast wave F^mc^ = E^j = const, and rewrite 
as 

= 0.3 (1 + z) 9u Jy, (82) 

which gives 

f° 0.3 / l3oba \i/3 [ (1 + i/< i/„ < 

Jy ^ D| bo“HzJ \ (£)410“)W4 /«£j,(i+^),/6(>/« v < v ,< v „. 

(83) 

Here fobs is in seconds and no is in cm”^; tjj is given by equation (56). Same results were 
derived previously (e.g. Sari et ah 1998; Panaitescu & Kumar 2000). 


6.2. Pair-dominated part 

The characteristic frequency Facc(tobs) is given by 
Farr = 4 X 10^"^ 
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We have fy ~ const in a broad range of frequencies around z/acc (see § 5.2.2) and therefore 

Lt - const (85) 

The observed flux is given by 
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d3 

m^cc = Y J?acc ^0 


^P, (87) 

and Face is given by equation (10). Finally, we substitute equation (76) for F(tobs) and find 
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Equation (88) shows that before tdec the observed flux evolves as F^(tobs) = const 

~l/2 

where fi, changes very slowly while increases as fobs if the postshock magnetic flux is 

1 /Q 

conserved (eq. 77) and reaches {Rdec/Ra.cc)c^ at fobs = tdec- After fdec, the observed flux 
decreases as with const = {Rdec/R s.cc)cb^] then 

/ , N -3/4-(9/8)a 

F^ocfV.oc ^ , (89) 

where a ~ 0.2 is the slope of spectral function f^y, and we get the approximate decay rate 
oc fobg®. A faster decay of is possible only if the magnetic held is destroyed in the 
blast wave and e^acc decreases with R. 

6.3. Radiative cooling of 

The above analytical calculation of F^ assumed that the pairs are slowly cooling. We 
now check this assumption. 


6.3.1. The cooling cutojf of the pair distribution 

The radiative cooling dehnes the cutoff of the distribution 7 c(m, R) (§ 3.3) and the 
corresponding cutoff frequency of the synchrotron spectrum z/c(m, R), 

~ eB 

z/c(m, R) = 0.2r- 7 c(^) R)^ (90) 

nieC 

where T is the current Lorentz factor of the blast wave and B = B{ni,R). The radiative 
cooling does not affect L^{R) if idjn, R) > ■ymi'cn, R) ed, m ruacc- 

Both synchrotron and inverse Compton losses can affect the value of 7 c at rriacc- They 
are proportional to the magnetic energy density f/s(macc,R) and soft radiation density 
Us{m^mR), respectively. 

As the blast wave expands from Race to some R < Rdec, its energy density U remains 
approximately constant while egacc grows as (R/Racc)^ (assuming that the magnetic flux is 
conserved, see eq. 77), so ^^(macc, R) ~ {R/Rb-ccYcbU. Us{macc, R) is dominated by synchro¬ 
tron radiation transported from the forward shock, which is in the fast-cooling regime. The 
density of synchrotron radiation at the current position of the shock is Us{rh, R) ~ {ceCBY^'^U 
(the rest of electron energy is taken away by inverse Compton scattering of synchrotron ra¬ 
diation, e.g. Sari & Esin 2001). A fraction ~ 0.2 of this radiation is transported to the 
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downstream region of the blast wave.® Then the density of soft radiation at m ~ rriacc is 


U^[m^,R) ^ 0.2 


( 91 ) 


This implies that Us{niacc, R) is approximately constant R < Rdec- 

So, we hnd that the relative contribntion of inverse Compton scattering to the cooling 
rate at m ~ rriacc is given by 


C (niacc) R^ 


Usjm^ccR) _ Q 2 =02 (^ ^ 

lls(macc,H) \ R J \^bJ 


( 92 ) 


where the last eqnality assnmes conservation of magnetic flux. Below we consider blast waves 
with C(mace,-Rdec) ~ 1 ; then the cutoff at m ~ mace is shaped by synchrotron losses rather 
than inverse Compton scattering of the forward-shock radiation. 

The cooling peaks at the deceleration radius and shapes the cutoff of the spectrum, 
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where R ~ i?acc is the radius where the pair-loaded mass m ~ mace was shocked and 7 ( 1 -1-/5) 
is the preacceleration factor at this radius. The last equalities in these expressions assume 
magnetic flux conservation. In particular, the factors m/mdec and [ 7(1 -|- /3)]“®/"^ appeared 
in the expression for udm, Rdec) because of the flux conservation; they would be absent if 
eB^m, R) = const. 

When the blast wave reaches the deceleration radius, its pressure starts to decrease, 
and the afterward evolution of the distribution function at m ~ mace is fully determined 
by adiabatic cooling — radiative cooling is slower and has a negligible effect. The adiabatic 


®The transport coefficient is evaluated assuming isotropic emission from the forward shock in the blast- 
wave frame. Its dependence on the downstream Lagrangian coordinate m is weak as long as m <C m. For 
example, the contact discontinuity (m = 0) receives ~ 0.15 of the forward-shock emission. 
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cooling of the blast wave expanded from i?dec to a current R is described by the factor 
A = (r(§ 3.3), and the resulting cutoff Lorentz factor is 


7c(m, R) 



1/2 

Rdec) 1 


R> R. 


dec* 


(95) 


The corresponding cutoff frequency evolves as Uc oc which gives 
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It can be expressed as a function of fobs using equation (76) for F and a prescription for the 

1 /2 

magnetic held evolution in the blast wave. If the magnetic hux is conserved, evolves as 
(eq- 77), and we get. 


Z/c(n7', tdec) V^dec/ V^dec/ 


^obs ^ ^dec • 


(97) 


If cb remained constant at all m and R^ one gets the cutoff frequency vR'n^Rohs) oc 
almost the same as with hux conservation. 


6.3.2. The slow-cooling condition 


A given shell m in a blast wave of radius R is in the slow cooling regime if 7 c(^^, R) > 
7 m(jTi, i?). Using equation (22), 


7m(m, R) = 




m 


^ 7(1 + /3) m, 

we hnd with the hux conservation assumption. 
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(99) 

Note that no and Fq cancel out from this expression. This ratio remains the same at i? > i?dec 
as both 7 m(^, 7?dec) and 7 c(ni, i?dec) evolve adiabatically and are reduced by a common factor 
A. The ehect of the 7 -ray precursor enters through the factor Z“^[ 7 (l+/ 5 )]“^/^, which reduces 
the ratio. 
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We are interested in m ~ rriacc because peaks near rriacc- Equation (99) shows that 
the slow-cooling assumption is valid if 
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( 101 ) 


7. Comparison with the reverse-shock model 


We focused in this paper on the forward shock of the GRB explosion. Early optical 
radiation is also expected from the reverse-shock in the ejecta if their energy is dominated 
by baryons rather than magnetic held (Meszaros & Rees 1993; Sari & Piran 1999). Explosions 
with the reverse and e^-loaded forward shocks should have two soft emission components, 
and it is instructive to compare them. 


Emission from the reverse shock peaks when it crosses the ejecta, which occurs at the 
deceleration radius. Then the explosion Lorentz factor P decreases and the reverse-shock 
emission decays. The decay rate can be calculated as follows. The Lorentz factor of electrons 
emitting at a hxed frequency z/obs = (1 + 2 :)zz is 

1/2 

OC (5ejf)-'/2^ (102) 

where Rej is the magnetic held in the ejecta at the current radius R. The number of electrons 
emitting at frequency u is 

N. ^ (103) 

where W is the total number of shock-accelerated electrons and is assumed. With 

decreasing P, W is reduced for two reasons: (1) ■ji, increases and (2) y^ decreases as the 
whole nonthermal spectrum is shifted to lower energies by adiabatic cooling. The adiabatic 
cooling factor is R oc where P = (4/3)r^ho/ie^pC^ is the blast-wave pressure, and hence 

W OC y-^’+^RP-^ OC (Rejfy^ (104) 
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The observed synchrotron hux is given by 
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If no destruction of magnetic flux takes place, -Bej evolves as 

Bej OC fiej-R oc R, 

where 7 is the adiabatic index of the ejecta material. Then, 

_pRS ^ ^^P+1)/27 +(p- 1)/4 p(p+l)/^+p ^(p+i)/2 


(106) 

(107) 


In the case of a uniform ambient medium, ho = const, we have after the deceleration radius 
f oc tob/* and R oc which yields 

pRS ,—3 (p+ 1)/87—p/4+1/8 f ^obs ’ h 5/3, p 2.5 /ins') 

7 = 4/3. p = 2.5 

A destruction of magnetic held in the ejecta could only steepen the decay. 

The mechanism of reverse-shock emission is similar to that of the afterglow: a shell 
of material (ejecta in the case of reverse shock and m^cc in the case of loaded forward 
shock) is heated at i? < i?dec with a low energy per particle, and after i?dec the shell cools 
down passively (adiabatically) producing a decaying hux of soft emission. 

The main difference between the two cases is the spectrum of emitting particles. The 
reverse-shock is thought to produce a power-law electron distribution with p = 2 — 3. By 
contrast, the effective spectrum of has p ^ 1. In the narrow shell Am ~ macci there are 
approximately equal numbers of at all energies up to the cooling cutoff, which gives the 
effective p about unity. The flat distribution of leads to the slow decay of their synchrotron 
emission that we found in § 6. The resulting light curve can, however, be made steeper if the 
ambient density decreases with radius or magnetic held is gradually destroyed downstream 
of the shock. Therefore, the main intrinsic difference between the reverse-shock and 
emission components is not the produced light curve but the slope of the distribution p. 

This difference can be observed directly by measuring the instantaneous synchrotron 
spectrum. The radiation is expected to have a small spectral index |q;| < 0.2 while the 
reverse-shock spectrum is much steeper: a = {p—l)/2 = 0.5 — 1 for p = 2 — 3. A measurement 
of the instantaneous spectrum in UV, optical, or IR at fobs ~ 100 s would provide a test of 
the current theoretical picture of the GRB explosion. Such a test can be done by Swift. 

Finally, we note one more difference between the reverse-shock and emission compo¬ 
nents. Both are cut off when < v, which happens at different times because the cooling 
frequencies Vc are different in the two cases. The ejecta magnetic held is likely stronger, 
Bej > R; then its cooling frequency is lower by the factor {Bej/B)~^ and hence the cutoff 
should occur sooner. 
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8. Discussion 


8.1. component of GRB afterglow 


An explosion blast wave is composed of swept-up layers (shells) of external medium that 
have been shocked in the forward shock front. In GRB explosions, these layers contain 
pairs injected into the external medium by the y-ray front. The layers shocked at small radii 
R < Rioad are dominated by the numerous nonthermal e^. At i? > i?ioad, these layers remain 
in the blast wave and form a thin shell adjacent to the contact discontinuity if there is 
no turbulent mixing in the blast. This shell has a low energy per particle and emits much 
softer synchrotron radiation compared to the outer swept-up material. The optical emission 
of the expanding blast wave can be strongly dominated by the shell even at large radii 
where its mass mioad is small compared to the total swept-up mass. If the shock magnetic 
parameter is below a critical value 10“^ — 10“^, the optical-emitting are in the slow- 
cooling regime and radiate their energy slowly, on a timescale longer than the deceleration 
time of the explosion. Their effect on the observed afterglow is described by adding a new 
emission component F^. 

The afterglow component is less sensitive to the model assumptions than the custom¬ 
ary pair-free afterglow. This is because the shell contains approximately equal numbers 
of particles with vastly different energies, which is a consequence of the steep evolution of 
the y-ray front and loading with radius. This fact allowed us to derive a simple formula 
for the observed spectral flux. 


V 


6 X 10“^ f2 


(1 + Z) ~l/2 „3/2 
n2 ^Ba.cc"'0 
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fu Jy. 


(109) 


It depends on the blast-wave Lorentz factor T that approximately corresponds to a given 
observer time fobs- In a broad range of u, the flux is dominated by a specihc mass shell 
m ~ rriacc whose shock radius was Race = 7 x cm. F^ depends on the magnetic 

held in this mass shell, Race which we parameterize using the usual equipartition parameter 
e_Bacc (Racc/S"^ divided by the energy density of the blast wave). Accurate calculations in §§ 4 
and 5 give the correction factor equation (109). All the details and uncertainties 

of the blast wave physics are absorbed by this factor (which is never much different from 
unity) and e^acc- 

To the hrst approximation, F^, decays after the deceleration time as T^ which is propor¬ 
tional to for an adiabatic blast wave in a uniform medium. The decay can, however, 
be faster if Race is gradually destroyed as the blast wave expands. 
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8.2. Uncertainties in the blast wave physics 

8.2.1. Magnetic field downstream of the shock 

The magnetic field behind the shock front and its downstream evolution are a major 
uncertainty of the current blast-wave models. In this paper we assumed, as is customary, 
that the magnetic parameter cb immediately behind the shock front remains constant as the 
shock expands, and followed the downstream evolution of eB{m,R) assuming conservation 
of the postshock magnetic flux. This assumption unambiguously determines the magnetic 
field in the shell, and we found that the shell is in the slow-cooling regime if the shock 
parameter cb is below a critical value 10“"^ — 10“^. If cb is above this value, the shell is 
in the fast-cooling regime and most of its energy is radiated at fobs < Uec- In this case, the 
radiation can be visible at fobs > Uec only in the infrared band; the time of the emission 
cutoff depends on frequency as 

The results are qualitatively similar if does not evolve in the postshock region ac¬ 
cording to the flux conservation but instead remains constant (which requires a gradual 
destruction of the magnetic flux with increasing distance behind the shock). However, the 
results would change if the magnetic flux is destroyed more quickly. The luminosity of the 
shell is proportional to its magnetic field and the field destruction would give a faster 
decay of Ffi. 

A real postshock field can be inhomogeneous on small scales (Medvedev & Loeb 1999) 
and one may need to include a distribution /(e^) in future afterglow models. This can 
change details. For instance, the cooling cutoff would be less pronounced because the low-e^ 
fraction would continue to radiate in the slow-cooling regime and give a tail of emission even 
when most of the postshock plasma is cooled and does not emit at the observed frequency. 


8.2.2. Turbulent mixing 

The e^-dominated material may not form a distinct thin shell near contact discontinuity 
if there are large-scale turbulent motions that mix up the postshock layers. Rayleigh-Taylor 
instability can drive such mixing like it does in supernova remnants. The mixing is unlikely 
to change the results of this paper. As long as we approximate the blast wave as a constant- 
pressure shell, the pressure in an e^-dominated gas element does not depend on its position 
within the blast wave. Therefore, its adiabatic and radiative cooling is the same as in the 
absence of turbulent mixing, and conservation of magnetic flux gives the same magnetic field 
in the element. 
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8.2.3. Mechanism of electron acceleration 

A significant uncertainty in the afterglow physics is the mechanism of electron accel¬ 
eration. A preshock magnetic held is compressed in the relativistic shock and becomes 
transverse, and the helds generated by Weibel instability are also transverse — parallel to 
the shock plane (Medvedev & Loeb 1999). Standard diffusive acceleration is unlikely to 
be efficient under such conditions. It requires the electron to cross the shock front many 
times, which can hardly happen: the upstream diffusion across the transverse held is slow 
because the electron gyroradius is smaller than the front thickness (proton gyroradius), and 
the electron will be advected downstream with the how velocity c/3 with respect to the front 
before it gets a chance to dihuse back to the upstream region. 

An alternative mechanism is stochastic acceleration by turbulence downstream of the 
shock. It may not be well described as impulsive acceleration and may keep electrons ener¬ 
getic even in the presence of signihcant radiative or adiabatic losses. This could change the 
theoretical afterglow light curve. 

However, the main signature of emission — white spectrum — will likely persist 
because the blast wave will still have a steep variation of density with the Lagrangian 
coordinate m, which invariably leads to the broad distribution with p ^ 1. This special 
feature ultimately comes from the y-ray front evolution with radius and is not related to the 
mechanism of acceleration. 


8 . 2 . 4 . Electron distribution function 

Customary afterglow calculations assume an idealized distribution function of the post¬ 
shock electrons: all electrons reside in a power-law component that starts at 7 ^ and ends 
at 7 c. In reality, it is possible that only a small fraction of electrons Ce are accelerated in a 
shock wave and the rest of them form a quasi-Maxwellian distribution; this is observed to 
be the case for collisionless shocks in the solar system. The expected energy of the accel¬ 
erated electron population Cg is then typically ~ 1% of the total plasma energy (which is 
dominated by the hot ions). Similar electron acceleration was inferred for supernova shocks. 
By contrast, Cg inferred from the existing fits of GRB afterglows is Cg ~ 0.1 (e.g. Panaitescu 
& Kumar 2002). 

We point out that the high Cg may have been inferred because the fits assume Ce = 1- 
Equally good hts may be obtained with a more reasonable ~ 0.1 and Cg ~ 0.01. This 
can be understood by looking at how the parameters enter the emission model. Besides 
the spectral slope p, the emitting electrons are described by two parameters: the number 
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of accelerated particles Ne = Ce-^t (where Nt is the total number of swept-up particles) 
and the minimum nonthermal Lorentz factor 7 ^- Note that is comparable with the 

mean energy per electron, which is proportional to Ce/Ce- Relaxing the assumption Q = 1, 
one can get the same 7 ^ by decreasing oc Q, and then the same observed emission may 
be explained with 10 times lower, and physically more plausible, values of C,e and Cg. The 
number of emitting particles Ne will not be changed if Nt is increased by the factor 
Thus, the reduction of Ce implies a higher ambient density and the circumburst density may 
have been systematically underestimated in the afterglow hts by one order of magnitude. 

Different assumptions concerning the shape of electron distribution function can lead to 
different afterglow radiation. However, the afterglow component is almost insensitive to 
such assumptions. For illustrative purpose, consider an extreme case. Suppose there is no 
power-law acceleration at the shock front, and the shock produces a narrow distribution 
peaking at 7 e(R) ~ ee{mp/me)T{^Z)~^. The factor 7 Z is determined by the 7 -ray transfer 
through the ambient medium and has a robust steep dependence on R near Race (§ 2 ). 
Therefore, 7 e will depend steeply on R, and the resulting distribution in the shell Am ~ 
mace swept-up near Race will be broad and flat. The same formula (eq. 109) will describe 
the radiation with a slightly different numerical factor ~ 1 . 


8.3. Prospects of detection of emission 

The radiation is predicted to have a special feature that can be tested with upcoming 
observations: its spectral index a is close to zero (|q;| ^ 0.2). This is a signiheant difference 
from the reverse-shock model which predicts a steeper spectrum, a = (p — l)/ 2 , where 
p = 2 — 3 is a putative slope of the electron distribution. Although the early optical emission 
has already been caught in a few bursts, no spectral data are presently available. Swift can 
provide the valuable spectral information. 

The dominance of the early UV/optical/IR afterglow can result in a characteristic 
two-peak shape of the light curve (see an example in Fig. 3). The emission component 
begins to decay at the deceleration time (eq. 74), which can be before Swift detects the 
afterglow. However, its decay is relatively slow {Fjf oc F^ oc t^^^* in the first approximation) 
and its tail is observable on timescales of minutes until the pair-free component takes over. 
The radiation should be visible for a longer time in the infrared band. 

To date, early optical emission (fobs < 1000 s) has been detected in four bursts: 
GRB 990123, GRB 021004, GRB 021211, and GRB 030418. In only one of them, GRB 990123, 
the peak of the optical flash was caught (Akerlof et al. 1999). This peak overlapped with the 
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prompt MeV burst, and hence the model developed in the present paper does not apply to 
GRB 990123: the must be Compton cooled by the MeV photons (keV in the fluid frame) 
and most of the energy is likely emitted in the GeV-TeV band (Beloborodov 2005). The 
strong optical radiation is expected in bursts where the MeV radiation front completely 
overtakes the blast wave by the time it reaches i?acc as discussed in § 1. 


8.4. Neutron front 

We studied in this paper the pair-loading effects on the forward shock. If the ejecta 
contains baryons, i.e. is not a pure electro-magnetic outflow (Poynting flux), a significant 
fraction of the baryons must be free neutrons (Derishev, Kocharovsky, & Kocharovsky 1999; 
Beloborodov 2003a). The neutron ejecta get completely decoupled and coast freely by the 
beginning of afterglow emission with a Lorentz factor T^ ~ Tej. They gradually /5-decay, 
however, some neutrons survive till radii R ~ 10^^ cm when the blast wave may have already 
decelerated, overtake the blast wave, and decay ahead of it depositing significant momentum 
and energy into the ambient medium (Beloborodov 2003b). Thus, GRB explosions are likely 
to develop leading neutron fronts that change the mechanism of the blast wave. 

The neutron front may emerge either after the e^-loading, at i? > Rioad, or at smaller 
radii, depending on the presence of fast neutrons with Tn > T. The impact of a fast neutron 
front on the early afterglow was recently studied by Fan, Zhang, & Wei (2004), and an 
alternative scenario with slow neutrons was proposed by Peng, Konigl, & Granot (2004). 
Accurate afterglow calculation that includes both the neutron decay and loading is a 
challenging problem which can be solved in future. We expect that the main signature of 
e^-loading — soft emission with a broad flat spectrum — be present in neutron-fed afterglows 
as well. 
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Appendix: Calculation of integral 


We here calculate the integral neglecting the cooling cutoff of the distribution 
function (i.e. assuming z>c > Then, 
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. 7(1 + /3). 

We have written 5L^^^/5m in this form to separate quantities that do not depend on the 
Lagrangian coordinate (functions of the current radius R only) and are constant in the 
integral. In a similar way, we rewrite the expression for Pmijn, R) (eq. 65), 
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Since lymifn) = R) is a monotonic function of the Lagrangian coordinate 0 < m < 

mi, we can change the integration variable to z/^. 
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where z>i is the peak frequency of synchrotron emission at mi, 
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The partial derivative of can be written as 
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varies with m much faster than e]/^rno ^'^^(1 + /3) (in the latter, only T could vary 
signihcantly, but even that does not happen for explosions with i?dec > -Race)- Therefore, the 
second term in equation (115) is dominant and the first term can be neglected. Then, 
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Z /and 7 are given as fnnctions of ^ in equations (4) and (5), and we find 

f dVrn \ 10^(H/^acc Ul > TTlacc, 

\ dm J m k \l7/2 m < mace- 


(117) 


Using equations (110) and (111) we get 
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Z needs to be expressed as a function of Vm- From equations (4) and (5), 

_ r 1 Z//ie < eV2 ^ 74, 

“ 1 {2Zle^piey/^ Z/Pe > eV2, 

which we substitute into equation (112) and then express Z from that equation, 
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The quantities no, e^, and F vary with m and therefore vary with z/^ when z/^ is chosen as 
the independent variable. However, their variation with z/^ is slow (because z/^ is a steep 
function of m). Therefore, with sufficient accuracy, we have Z oc at m{vrn) > mace and 
Z oc at m{Vm) < mace- 

Now the integral (118) can be calculated. We consider hrst R < R^cc and then R > Race- 


R — Race 


At radii R < Race the Lagrangian coordinate m < m < mace] then dlnum/dlnm = 
{17/k) and Z oc 

For z/ > z>i we have z/^ < ^ in the whole blast, and the integral (118) reads 
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where Zi = Z{mi) and we have used Z/Zi = (z/m/f'i)”^^^'^- We also used the fact that 
and m vary slowly with z/^ and took them as constants evaluated at z/^ = i^i where the 
integral peaks. 



50 


For h> < h>i the integral can be written as a sum of two integrals over Um < and Um > v. 
Both integrals peak at We denote m* = m(z/m = 

use ZjZ^ = {vm/I and get 
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We used again the slow variation of m and Ib with z/^ and replaced m by m* and by eg*. 


R 'Z. RqjCC 


At radii R > Ra,cc, the blast wave material has shells with mass coordinate 0 < m < rriacc 
and rriacc < m < fh. The integral is taken over 0 < m < mi = min{?7i, mioad}- 

For z/ > hi the integral is given by 
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Let us denote 

^acc(7?) = Z/jyj(macc) 7?). (124) 

We have Z oc z/m at z/^ < t'acc and Z oc z/m at z/^ > t'acc- It is convenient to calculate 
the integral as a sum of two integrals over < hacc and z/^ > 7’acc- The hrst integral peaks 
at Vyn = hacc and the second at z/^ = T’l- We then hnd 
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(125) 


where Zacc = ^(^acc) = 74/ie, esacc = eB(^2racc,-R), and eui = eB(mi,.R). As soon as 
the blast wave radius exceeds i?acc, we have almost immediately hacc -C hi (hi increases 
exponentially with R between i?acc and i?ioad)- Then the terms with hacc/hi can be neglected 
(these terms are needed only to match smoothly the formula with the results obtained at 
R < 7?acc)- 

Next, consider intermediate v in the range hacc < z/ < hi. Then the integral is calculated 
as a sum of three integrals over 0 < z/^ < T’acc, f'acc < and z/ < z/^ < hi. In the first 
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interval we use ZjZ^c = (t'm/^'acc) and in the other two intervals ZjZ^ = {vmlv) 
where asterisk refers to the point = v. We thus get, 
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Finally, at low z/ < hacc, the integral is calculated as a sum of three integrals over 
0<h'^<u,u < 12 ^ < z>acc, and z/acc < < f'l- In the hrst two intervals we use 

Z/Z^ = ^ and in the last interval ZjZ^^. = {i^ml^a,cc)~^^‘^■, 
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Z* and m* 


The above formulae for should be made explicit by expressing m*, .^*, and Z* in 
terms of z/, R. 


Equation (120) gives Z(z/m), and we get Z = Z* at z/^ = n, 
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All quantities with asterisk are taken at i? = i?* (z/^ = n). Since i?* is very close to R^cc 
(see eq. 133 below), it is sufficient to use the first approximation i?* = i?acc in the equation 
for Z*. Then we get, 
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and the subscript “acc” marks that the quantity is taken at m = rriacc- 
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The relation between ^ and Z is given by equation (4), and we hnd 
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(here <h = 8 x 10^ corresponds to = 3.^acc)- Using the relation between R and ^ (eq. 3), 
we get 
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where we took $ + — 1 ~ 2$. The corresponding m* is 
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Final result 


1 /2 

The terms with are non-negligible when u ~ z/acc and m* ~ niacc, therefore one 
can replace by and simplify the derived expressions. The hnal result is as follows. 
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In this paper, we are interested in i? > i?acc and v < Vi. Then 

= 30 Q^m^Z^, 

and the numerical factor fi, dehned in equation (67) is given by 


This factor is shown in Figure 9 assuming = e^acc- 
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